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ABSTRACT 

We present the first results of a survey for high redshift, z ^ 6, quasars using izY multi-colour 
photometric observations from the Dark Energy Survey (DES). Here we report the discovery 
and spectroscopic confirmation of the zab, Yab = 20.2, 20.2 (M1450 = —26.5) quasar DES 
J0454—4448 with a redshift of z = 6.09±0.02 based on the onset of the Lyman-G forest and 
a HI near zone size of 4.11^ 2 proper Mpc. The quasar was selected as an i-band drop out 
with i—z = 2.46 and zab < 21.5 from an area of ^300 deg^. It is the brightest of our 43 
candidates and was identified for spectroscopic follow-up solely based on the DES i—z and 
z—Y colours. The quasar is detected by WISE and has WIab = 19.68. The discovery of 
one spectroscopically confirmed quasar with 5.7 < z < 6.5 and zab ^ 20.2 is consistent with 
recent determinations of the luminosity function at z ^ 6. DES when completed will have 
imaged ^5000 deg^ to Yab = 23.0 (5cr point source) and we expect to discover 50-100 new 
quasars with z > 6 including 3-10 with z > 7 dramatically increasing the numbers of quasars 
currently known that are suitable for detailed studies. 

Key words: dark ages, reionization, first stars — galaxies: active — galaxies: formation — 
galaxies: high redshift — quasars: individual: DES J0454-4448 


1 INTRODUCTION 

High redshift quasars act as important probes of the very early Uni¬ 
verse. They can be used to study the formation and evolution of 
the first supermassive black holes and to provide insight into the 
properties of the inter-galactic medium (IGM). The Lyo forest and 
associated absorption lines blueward of the Lyo peak in a high- 
redshift quasar spectrum give a direct measurement of the neutral 
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hydrogen in the IGM (Fan et al. 2006). Detailed studies of the ion¬ 
ized near zone in quasars can be used to investigate the ionization 
state and evolution of the neutral hydrogen fraction at z>6 (Bolton 
& Haehnelt 2007). Cosmologically distributed intervening quasar 
absorption lines allow us to map the distribution of gas in the IGM 
as well as the metal content. Identifying larger numbers of high 
redshift quasars would enable such studies to be conducted further 
back in cosmic time along independent lines of sight through the 
IGM. 

Spectra of quasars from z ~ 2 (Gunn & Peterson 1965) up to 
z ~ 6 Becker et al. 2007; Fan et al. 2006 indicate that the IGM 
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is highly ionised {nm/nii ^ 10“^) and therefore that reioniza¬ 
tion was complete at earlier epochs. The 2cr lower-bound for the 
Epoch of Reionization is estimated to be at z ~ 6 from Cosmic Mi¬ 
crowave Background (CMB) measurements (Planck Collaboration 
et al. 2015). The discovery of a large sample of new high redshift 
quasars at z > 6 would allow this important epoch to be studied in 
detail. Even the discovery of a few new objects can be valuable if 
they have extreme properties, such as being very bright ((Wu et al. 
2015)) or very distant (Mortlock et al. 2011). 

Previous surveys for high redshift {z ^ 6.0) quasars (e.g. Car- 
nail et al. 2015; Pan et al. 2006; Jiang et al. 2009; Mortlock et al. 
2012; Venemans et al. 2015, 2013; Willott et al. 2010) have already 
led to the discovery of ~ 50 quasars at z ^ 6.0. Most of these sur¬ 
veys have employ purely optical photometric datasets such as those 
provided by the Sloan Digital Sky Survey (SDSS) and the Canada 
Prance Hawaii Telescope Survey (CPHTLS). The new optical pho¬ 
tometry provided by the Dark Energy Survey (DPS) (The Dark En¬ 
ergy Survey Collaboration 2005) provides one notable advantage 
over most of these previous optical surveys from the point of view 
of high redshift quasar selection. The Dark Energy Camera (DE- 
Cam) CCD detectors are optimised to extend to the near infrared 
wavelengths of ^l/xm. The associated grizY filter-set therefore 
allows detection of quasars out to higher redshifts than was possi¬ 
ble with previous optical surveys. The DPS z—band and the red¬ 
der y—band coupled with the large area of the celestial sky not 
previously surveyed and the increased depth over previous surveys 
enables detection of Lya at 6.5 ^ ^ ^ 7.2. In previous optical sur¬ 
veys the sensitivity in the z band resulted in making it increasingly 
difficult to find quasars with increasing redshift. 

Here we present the first search for high-redshift quasars in 
the DPS Science Verification Data. DPS magnitudes, near infrared 
(NIR) VISTA magnitudes and WISE magnitudes are quoted on the 
AB system. The conversions from Vega to AB that have been used 
for the VISTA data are: Jab = Jvega + 0.937 and Ksab — 
Ksvega +1.839, thcsc arc taken from the Cambridge Astronomical 
Survey Unit’s website. ^ The conversions for the ALLWISE data 
are WIab — Wlvega + 2.699 and W2ab — W2vega + 3.339 
which are given in Jarrett et al. (2011) and in ihd ALLWISE explana¬ 
tory supplement. ^ When required a flat cosmology with Dmo = 0.3 
and Ho = 70.0 was used. 


2 DARK ENERGY SURVEY DATA 

The Dark Energy Survey is a 5000 deg^ optical survey in the south¬ 
ern celestial hemisphere being conducted using the Dark Energy 
Camera (DECam) (Plaugher et al. 2012). DECam is mounted on 
the Blanco 4-meter telescope at the Cerro Tololo Inter-American 
Observatory (CTIO) and the data are processed by the DPS data 
management system (Desai et al. 2012; Mohr et al. 2012). The 570 
Megapixel DECam has a field of view of 3 deg^ with 0.27” per 
pixel. To allow for a good detection efficiency of high redshift ob¬ 
jects, DECam CCDs are very sensitive to the red part of the spec¬ 
trum. These fully depleted, 250 gm thick detectors have been de¬ 
veloped at the Lawrence Berkeley National Laboratory to be used 

^ http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/filter-set 
^ The ALLWISE explanatory supplement, 

http://wise2.ipac.caltech.edu/docs/release/allwise/expsup/sec5_3e.html, 
directs the reader to the WISE All-Sky ex¬ 
planatory supplement for the conversions; 

http://wise2.ipac.caltech.edU/docs/release/allsky/expsup/sec4_4h.html#summary. 


for DPS. The quantum efficiency of these devices in the z band is 
above 50%, almost an order of magnitude higher than traditional 
thinned devices. This also allows Y band observations to be taken 
within the same survey reducing the need for follow-up photome¬ 
try to confirm high redshift quasar candidates. The combination of 
depth, red sensitivity and area provided by DPS, makes it an ideal 
experiment for searching for luminous quasars at high redshift in 
the southern celestial hemisphere. 

The DPS survey area also overlaps the near infrared VISTA 
Hemisphere Survey (VHS) (McMahon et al. 2013) at ~l-2/im 
and the all-sky Wide Infrared Survey Explorer {WISE) data at 3.4, 
4.6, 12 and 22gm. These two surveys therefore provide infrared 
photometry for high-redshift quasar candidates selected in DPS 
(Banerji et al. 2015). 

The DECam obtained its first light images in September 2012 
and this was followed by a period of Commissioning and Science 
Verification observations (SV hereafter) lasting between October 
2012 and Pebruary 2013. In this work we utilise the Science Veri¬ 
fication Pirst Annual Data Release (SVAl) by DPS Data Manage¬ 
ment to conduct our high-redshift quasar search. The SVAl ob¬ 
servations were designed to reach the full DPS survey depth ex¬ 
pected after the nominal 5 year survey period. However, the depth 
across the SVAl area is non-uniform. Pigure 1 shows the cumu¬ 
lative area in the DPS SVAl region down to a given lOcr depth 
in the z-band derived from mangle (Swanson et al. 2008) based 
polygon depth maps. These lOcr depth maps are calculated for a 
2” diameter apertures. These values will depend on the seeing in 
each polygon and have been converted to PSP magnitudes using 
the median difference for all stellar objects (as defined in section 
X) The median of these differences (Zpsf — ^ 2 ”) was found to be 
-0.28 with a ctmad = 0.066. This offset was applied to the mag¬ 
nitude limit calculated as shown in figure 1. In this paper PSP and 
model magnitudes are used. The PSP magnitude is the point spread 
function magnitude and the model magnitude is the magnitude once 
an galaxy model has been fitted to pick up the fiux in the extended 
wings of the object. Both of these parameters are produced by the 
DES data management pipeline. As detailed later, a fiux-limited 
sample down to a depth of Zpsi < 21.5 is used for the high-redshift 
quasar search presented here, which corresponds to a total survey 
area of 291 deg^. 


3 QUASAR CANDIDATE SELECTION 

Whilst low redshift (z < 2) quasars have characteristically blue 
optical colours (from ultra violet excess) (Schmidt & Green 1983) 
at ^ > 4 quasars start to become red in the shorter wavelength opti¬ 
cal bands (Hook & McMahon 1998; Storrie-Lombardi et al. 2001). 
The observer frame optical colours of these high redshift quasars 
are characterised by absorption due to neutral hydrogen in the in- 
tergalactic medium in the rest frame UV shortward of the Lyman-a 
(Arest = 1216A) emission line. At z = 6 the Lya line is redshifted 
to 8512A which leads to red {i — z) colours for quasars at these 
redshifts. These high-redshift quasars are therefore selected in op¬ 
tical surveys as z—band drop outs (Pan et al. 2000). Below we detail 
our selection method for isolating high-redshift quasar candidates 
from the DES imaging data. PSP magnitudes are used throughout 
(unless otherwise stated) as they provide the best measure of flux 
for unresolved point-like sources. 

The primary colour selection used is an (i — z) colour crite¬ 
rion, which is based on the selection used by Pan et al. (2001). The 
Pan et al. (2001) SDSS selection of {i — z) > 2.2 was mapped 
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Figure 1. Cumulative area versus lOcr z-band depth in a 1 " diameter aper¬ 
ture and PSF magnitudes for the DBS SVAl data. The aperture magnitudes 
were converted to PSF magnitudes using the median offset between PSF 
magnitude and aperture magnitude for point sources from the whole SVAl 
dataset. Our magnitude limit of < 21.5 is shown as the vertical line. 
The aperture magnitude limits were taken from the DBS Mangle (Swanson 
et al. 2008) products. 


14 



Figure 2. The black points are all the objects on one tile from the survey 
and the red points show the quasar candidates after the final stage of the 
selection in section 3. A clear locus can be seen at 0.0 on the horizontal 
axis. The blue points are known quasars taken from the catalogue given in 
Veron-Cetty & Veron (2010) with 2 ; > 2.0. Below this redshift host galaxy 
contamination can occur. The median of the Veron 2010 sample is 0.015 
and the ctmad = 0.052. The solid vertical line shows the cut used to select 
point sources. The dashed horizontal line shows the magnitude limit used. 


onto the DBS filter system as follows. First we selected point-like 
objects in the area of overlap between the DBS and SDSS surveys, 
as those having Zpsf — ^modei < 0.145. This point source selec¬ 
tion is illustrated in figure 2 and is detailed below. The DBS and 
SDSS sources were matched using a radius of 2.16" (0.006°). A 
least squares linear fit was carried out on a running median of the 
colours. A window size of 50 was used for the median and the data 
was sorted by (i — ^)sdss and then sparsely sampled along the 
(i — z)sDSS direction with a stride of 10 to determine the colour- 
term between the two surveys: 

(i — z)des = 0.72(i — z)sDSS + 0.111 (1) 

Based on this and the Fan et al. (2001) selection criterion we derive 
an equivalent colour criterion of (i — 2:)des > 1.694. 

We now list each of the steps in our high-redshift quasar se¬ 
lection process, which are also summarised in Table 1. 

(i) High-redshift quasar and other rare object searches are of¬ 
ten contaminated by a large number of artefacts. These spurious 
sources often dominate outlier regions of colour space that are also 
populated by smaller numbers of real objects. In fact it is safe to as¬ 
sume that most outliers in an unfiltered catalogue are junk. The first 
step therefore involved artefact removal at the catalogue level. Ob¬ 
jects that were flagged by the DBS data management SVAl pipeline 
as having bright neighbours that were close enough to bias photom¬ 
etry, being originally blended with another object or having at least 
one saturated pixel were provisionally labelled as junk. If they were 
also flagged as having at least one low-weight pixel in the measure¬ 
ment image within the isophotal footprint or in the filtered detection 
image then this was confirmed and they were removed. A combina¬ 
tion of the weight flags and the image flags was used to ensure that 
faint objects undetected in some of the bands were retained while 
removing bright saturated objects and bleed trails. The flags cut re¬ 
sulted in 36,959,467 objects remaining of the original 40,129,963 
in the DBS SVAl release. Examples of the artefacts removed in this 
step are shown in Appendix Al. 

(ii) A flux limited sample was then created by limiting the sam¬ 


ple to sources with: 

Zpsf ^ 21.5 and cr^ < 0.1 (2) 

(iii) Objects were required to satisfy a point source criteria of: 

Zmodel ^ 0.145 (3) 

which is the star-galaxy classifier defined in the SDSS photomet¬ 
ric survey (Stoughton et al. (2002)) and SDSS quasar selection 
(Richards et al. (2002)). After this step there were 4,790,106 ob¬ 
jects remaining. This star-galaxy separation criterion is shown in 
figure 2. The stellar locus is clearly situated at Zpsi — ^modei ~ 0. 
We can see that the star-galaxy separation starts to become less reli¬ 
able at fainter magnitudes but the stellar locus is still clearly visible 
down to z < 21.5, which is the flux limit imposed on our sample. 

(iv) The colour criterion derived above: 

(i — 2:)des > 1.694 (4) 

was then applied leaving 23,517 objects. 

(v) As high redshift (z > 6) quasars should not have detectable 
flux in their rest frame shortward of the hydrogen limit limit at 
Arest = 912A corresponding to Aobs ~ 6500A they should be 
undetected in the g or r bands. Visual inspection of a sample of 
candidates indicated that some had significant detectable flux in g 
and r. We therefore applied the following criterion to the g and r 
band measurements 

^psf > 23.0 and rpsf > 23.0 (5) 

which were chosen as this ruled out almost all of the objects that 
were reliably detected without removing undetected sources in re¬ 
gions of particularly deep photometry. 

(vi) There were a few objects that were very faint in the g or 
r band but were in a very deep part of the survey and still had a 
reliable detection in these bluer bands. To remove these detected 
objects an error threshold of Ug and cTr >0.1 was used. i.e. Objects 
with cFg and Gr ~ 0.2 (5cr detections) are retained since there is a 
finite probability that they have zero flux. 
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(vii) Cool dwarf stars, such as L and T-dwarfs, in our own 
Milky Way, are the main astrophysical contaminants in high red- 
shift quasar searches. These were removed using a colour cut of 

2^psf ^ 0-5 (6) 

The colour selection can be seen in Figure 3 which shows the 
cool dwarf locus clearly distinct from the selection box in terms of 
the 2 : — y colour. We note that a blue z — Y cut will also potentially 
remove higher redshift quasars at z > 6.5, at which point the Lya 
peak begins to move into the DBS Y -band, making the 2 ; —y colour 
redder. In the future, we will use NIR data, particularly in the J- 
band from VHS, to push our quasar search to higher redshifts. This 
selection will cover 2 : — y dropouts, rather than i — z, and will be 
discussed further in future work. 

(viii) At this stage many of the objects left were bright artefacts 
in the z band which were not present in any of the other bands. 
To limit the numbers of these that was present in the sample two 
Y band magnitude cuts were used. First objects with bad Y band 
photometry (Ypgf = 99) were removed and then those with Ypgf < 
23.0. These cuts ensured that objects were detected in the Y band 
as well reducing their likelihood of being a z band artefact. 

(ix) Following these cuts, many of the objects that remained had 
PSF magnitudes of 99 in the r-band. They corresponded to two 
classes of objects: either very bright objects or those at the r-band 
detection limit. The bright objects were removed using list driven 
forced photometry in the r-band. For each object we identified the 
brightest pixel in a box of 3x3 pixels centred on the position of the 
object and required that it was no more than 10 analog to digital 
units greater than the median value in a 30" box around the object. 
This corresponds roughly to three times the median ctmad (3.06) 
above the median of the median values (0.147) of the images in 
the data release. A conservative threshold of 10 was used to filter 
these bright objects to ensure that any faint objects which may cor¬ 
respond to high-redshift quasars, were not removed. 

(x) Most of the objects that remained at this stage were cosmic 
ray detections in the z-band. These were removed based on the 
differences between adjacent pixels (as a proxy for the local pixel 
gradients) as a cosmic ray would be more abrupt in terms of its 
pixel gradient than a real object. In future DBS data releases, a new 
method of cosmic ray reduction will be included in the pipeline 
processing. 

(xi) Following these steps there were few enough candidates to 
cut down by systematic visual inspection. The visual inspection 
stage removed 34 spurious sources such as cosmic rays, satellite 
trails and diffraction spikes. At the end of this stage we have 43 
candidate high-redshift quasars. 

(xii) We ranked these 43 final candidates based on an error 
weighted photometric selection metric for follow-up spectroscopy. 
The selection metrics 0iz and 0 zy combine to give 0izY which es¬ 
sentially constitutes a error-weighted measure of how far into the 
colour-selection box, shown in Figure 3, each candidate lies. 


0iz 

0zY 

0izY 


(^psf -^^psf) 1.694) 

V^^err “f ^err 


0.5 ('S^psf ^^sf) 

.7^2 I y2 
V ^err i err 



4>zY 

|0zy| 


/>z\ 


(7) 

( 8 ) 
(9) 


For those sources undetected in the i-band, the i-band magnitude 
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Figure 3. The black points represent all objects that have 2 ;psf — 2;modei < 
0.145 and are detected in the i, z and Y bands from one tile of the survey. 
The blue points are the candidate objects that satisfy the first 11 steps of the 
selection criteria. The objects that pass the final step are circled while the 
others are faded. The blue star shows the location of the confirmed quasar. 
The grey line shows a simulated quasar track starting from z = 5 and 
incrementing by 0.1 to 2 : = 6.4. Objects that have i = 99 (are undetected 
in the i band) pass the selection criteria and are here plotted with arrows to 
indicate that they are lower limits in i-z. The objects outside the selection 
box are limits in the i band where the survey is particularly shallow. Many 
of the objects are situated at the lower right of the selection box which 
suggests that they are cool stars scattered in by random photometric errors. 
This is supported by the large errors that many of them have which overlap 
the selection box edges. The red, orange and brown points are known brown 
dwarfs found with UKIDSS and converted using the spectra and the DBS 
band passes. 


limit at the source position was used instead in the computation of 
0iz and the error set to 0.2. 

This selection led to prioritising four of the 43 candidates as 
highly viable high-redshift quasar candidates with 0izY > 4 from 
the 291 deg^ region. The distribution of 0 values can be seen in Fig¬ 
ure 4 where the four candidates have been circled. One of these can¬ 
didates clearly has much larger values of 0 than the other three and 
this object (DBS J0454—4448) with zab = 20.2 was therefore se¬ 
lected for follow-up spectroscopy as our most promising first high- 
redshift quasar candidate from the DBS data. DBS J0454—4448 is 
shown across the DBS and VISTA bands in figure 5. In these cutout 
images the absence of flux in the g and r bands can be clearly seen 
along with the increase in brightness from the i to z bands. 


4 DES J0454-4448 

The photometric data for DBS J0454—4448 is summarised in Ta¬ 
ble 2. This includes the DBS photometry (g, r, i, z and Y) as well as 
infrared magnitudes from the VISTA Hemisphere Survey (J and K) 
and the AllWISE data release of the WISE All-Sky Survey (Wright 
et al. 2010) (W1 and W2). Blain et al. (2013) have studied the in¬ 
frared WISE colours of known high-redshift quasars finding that 
they are typically bluer than the WIab - W2ab >0.16 selection 
cut used to select lower redshift AGN such as in Stern et al. (2012). 
The colour (WIab - W2ab = 0.05) of this object supports this 
and agrees with the locus of known high redshift quasars matched 
to the WISB data. The colour also agrees with the predicted track in 
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Figure 5. Cutouts of our spectroscopically confirmed quasar, DESJ0454-4448 in each of the DES (g, r, i, z and Y) and VHS (J and K) wavebands. The boxes 
are 30” across and the quasar is the central object. 



Figure 4. The (f)izY selection metric. Here the x axis is the (piz defined in 
equation 7 and the y axis is 4>zY • The four objects with (j)[zY > 4 are circled 
and the object that was confirmed as a quasar has the highest value of both 
^iz and ^zY the figure. 


Table 1. Summary of the steps in the high-redshift quasar selection process 


Step 

Description 

Number 

Removed 

Number 

Remaining 


Number of objects in database 


40,129,963 

1 

Flag criteria 

3,170,496 

36,959,467 

2 

Zpsf ^ 21.5 and fJz <0.1 

29,760,121 

7,199,346 

3 

Zpsf '2^model ^ 0.145 

2,409,240 

4,790,106 

4 

ipsf " Zpsf ^ 1.694 

4,766,589 

23,517 

5 

gpsf ^iid rpsf > 23.0 

12,985 

10,532 

6 

ag and ar >0.1 

371 

10,161 

7 

Zpsf " Ypsf ^ 9.5 

92 

10,069 

8a 

Ypsf not equal to 99 

8725 

1344 

8b 

Ypsf < 23.0 

506 

838 

9 

Forced photometry in r band 

372 

466 

10 

Cosmic ray removal 

489 

77 

11 

Visual Inspection 

34 

43 

12 

0izY > 4 

39 

4 


Assef et al. (2010) as shown in Figure 7. As the resolution of WISE 
(6.1” and 6.4” in W1 and W2 respectively) is quite large compared 
to the separation of the object and its nearest neighbour (4.89”) it 
is debatable if the magnitude corresponds solely to the quasar. 

Spectroscopic follow-up of DESJ0454—4448 was conducted 
using the Magellan Echellete (MagE) spectrograph (Marshall et al. 
2008) on the 6.5m Clay telescope at the Las Campanas Observa¬ 
tory in Chile. Two 10 minute exposures were taken on the night of 


Rest Frame Wavelength / Angstrom 

1180 1190 1200 1210 1220 1230 



Figure 6. In this plot the ratio of the continuum flux and the spectrum is 
shown. The red vertical line shows the position of the Lyman alpha peak. 
The dashed blue line is at a flux ratio of 0.1 and the dotted blue line show 
the zero level. The nearzone has a size of proper Mpc. 

December 2nd 2013 in variable conditions. The data were reduced 
and the one-dimensional spectrum combined using a custom set of 
IDL routines. The discovery spectrum for the quasar can be seen in 
Fig. 8. 

Measuring the quasar redshift from emission lines in the 
present spectrum is difficult due to the strong attenuation of Lya 
(Arest=1215.67A) and the modest quality and limited wavelength 
range of the data which does not cover the wavelenght range to 
detect CIV(1549) or SiIV-fOIV](1400) which at z ^ 6.00 have ob¬ 
served wavelenghthat A ^ 9800 A. We therefore rely on the appar¬ 
ent onset of hya forest absorption to estimate the redshift. This 
approach has been previously used for weak lined quasars such as 
BL Lac objects (Danforth et al. 2010). 

A series of absorption lines appears over 8530 < A < 8630 A, 
which we interpret as the expected Lya absorption in the quasar’s 
near zone. The reddest of these lines falls at 8626 A, which corre¬ 
sponds to a Lya redshift of z = 6.096. We tentatively detect the 
quasars Si II 1260 A emission line with a peak wavelength near 
8920 A, though it is contaminated by skyline subraction residuals, 
corresponding to a redshift of 6.077. Combining these estimates 
gives a mean redshift and standard deviation: z = 6.09 ±0.01. To 
refiect the systematic uncertainty in the redshift we follow Ean et al. 
(2006) and adopt an total uncertainty of 0.02 to take account of the 
observed offsets in redshifts between low and high ionization lines 
in quasars. 

The quasar continuum is strongly attenuated at A ^ 8530 A. 
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Table 2. Properties of DBS J0454-4448. The g band magnitude is given as 
a lOcr magnitude limit for a 2” aperture. 


DES J0454-4448 

DES Tilename 

DES0453-4457 

RA (J2000) 

73.50745 (04^54"^01.79^) 

DEC (J2000) 

-44.80864 (-44°48’31.1”) 

Redshift 

6.09 ± 0.02 

g 

> 23.41 

r 

25.30 ± 0.25 [AB] 

i 

22.66 ± 0.05 [AB] 

z 

20.20 ± 0.01 [AB] 

Y 

20.24 ± 0.04 [AB] 

J 

20.24 ± 0.07 [AB] 

Ks 

20.12 ±0.17 [AB] 

W1 

19.68 ± 0.08 [AB] 

W2 

19.62 ±0.14 [AB] 

Mi450 

-26.48 ± 0.03 

^NZ 

4.1 Ai Mpc 

Z,corrected 

4.8lJ;4 Mpc 


Following the method outlined in Fan et al. (2006), we take the 
edge of the proximity zone to be the point where the ratio between 
the continuum flux and the spectrum first falls below 0.1 blueward 
of the Lya peak. Here the intrinsic continuum is roughly estimated 
using a composite from low-redshift SDSS spectrum (as used in 
Mortlock et al. (2011)) normalized to the flux between 1268 and 
1295 A. The nearzone size for an edge at 8530 A and the quasar 
redshift estimated above is proper Mpc, the error in this 

measurement is dominated by the error in the redshift. In figure 6 
the normalised flux and the extent of the nearzone are shown. The 
corrected nearzone size of this object is found to be consistent with 
the measured relationship in the Carilli et al. (2010) study. Due to 
the low signal to noise of the data no constraints can be put on the 
optical depth at this redshift. 


5 DISCUSSION 

The number density estimates based on the Willott et al. (2010) 
luminosity function suggest that there will be between 50 and 100 
new quasars with z > 6.0 across the ~5000 deg^ of the full DBS 
area down ioY ab = 21.5. The predictions also suggest that there 
will be 3 to 10 quasars with z > 6.5 and Yab < 21.5 found in the 
complete survey footprint. 

Estimates from the Willott et al. (2010) and McGreer et al. 
(2013) luminosity functions suggest that there are 5 to 10 quasars 
with z < 6.5 in the SVAl area. Our selection method gives 43 
other candidate objects and suggests four as higher priority than 
the others. This final selection step is biased towards bright objects, 
due to them having smaller errors than fainter candidates, and is 
not designed to And a complete sample. To And all the quasars in 
this area with a high purity would require full Bayesian analysis 
of the quasar population and the brown dwarfs. The exact number 
of quasars predicted in our survey is also sensitive to the choice of 
quasar template. Different predictions for the hy—a forest in these 
quasar template spectra can change the exact redshift range probed 
by our colour selection criteria. This then heavily influences the 
predicted number of objects. 



Figure 7. The figure shows the evolution in W1-W2 colour with redshift 
of known quasars. The grey points are known low-redshift quasars from 
SDSS DRV with the red line being a median smooth of these data. This 
allows comparison with the predicted track, from Assef et al. (2010), which 
is shown by the black line. The blue points are a sub sample of the known 
high redshift quasars matched to WISE, these objects were matched to the 
ALLWISE data and then visually inspected to remove blended sources in 
WISE. The large blue star is DES J0545-4448. 


6 CONCLUSIONS 

We have presented the first search for high-redshift quasars in the 
Dark Energy Survey using 291 deg^ of optical imaging data from 
the DES SV observations. Our search algorithms are demonstrated 
to be effective in removing image artefacts and spurious sources 
such as cosmic rays and satellite trails that can populate the same 
regions of colour space as rare objects in the data. Based on the 
DES photometric data alone, we were able to reduce the number 
of candidates in this region to a size suitable for visual inspection. 
We also did not require confirmation photometric follow-up as has 
been the case in previous high-redshift (z > 6) quasar searches. 
A photometric error-weighted colour selection metric was intro¬ 
duced as a method of ranking candidates according to the prob¬ 
ability of being high-redshift quasars. The highest ranked candi¬ 
date, DESJ0454—4448 was spectroscopically confirmed to be a 
high redshift quasar at z = 6.10. Three more candidates were also 
identified as having colours consistent with high-redshift quasars 
and will be further investigated as part of the larger sample of high- 
redshift quasar candidates from the DES Year 1 observations. 

DESJ0454—4448 is also detected at infrared wavelengths in 
the VISTA Hemisphere Survey and the WISE All-Sky Survey and 
its infrared colours are entirely consistent with those expected 
for high redshift quasars. Using current estimates of the z ~ 6 
quasar luminosity function, we have shown that the first z ~ 6 
quasar is expected in the 291 deg^ region used in this work, at 
a z-band magnitude of 20.2. The measured z-band magnitude of 
DESJ0454—4448 is 20.2, entirely consistent with these predictions. 

We conclude that our search algorithms for high-redshift 
quasars are successful in isolating quasars at z > 6 from DES pho¬ 
tometry alone. The analysis method presented here should be easily 
scalable to the full 5000 deg^ area of DES and is therefore expected 
to result in substantial numbers of new high-redshift quasars in the 
southern hemisphere over the next few years. 
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Figure 8. The discovery spectrum of DBS J0454-4448 at z = 6.09. It shows the characteristic flux deflcit blueward of Lyman alpha indicative of a high redshift 
quasar. The red vertical lines over the spectrum show the positions of the Ly-/3 and the Ly-a emission line redshifted to z = 6.09 as well as the errors on the 
position of Ly-a from the error in the redshift. The middle plot shows the error spectra and the bottom plot indicates the efficiencies of the DBS Alters as a 
function of wavelength. Both the object spectra and associated error spectra have been smoothed to the same degree using a gaussian Alter. 
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Figure Al. All of the objects shown here are flagged with a three in the DES 
flags. It can be seen that all of these objects have very different problems 
with them and that the first cutout shows an object that we would want to 
keep in if it looked like that in the g or r band. These cutouts are 30” by 
30” and all taken from the z band images. The flagging was found to be 
similar across all the bands and objects shown here have many analogies in 
the other bands. 



Figure A2. In these cutouts an object is shown with ^fpgf = 23.16 ± 0.03 
and Fpsf = 23.02 ± 0.03. The source is obviously detected in the g and r 
band despite having a magnitude below our limit. 


2 The object was originally blended with another one 
4 At least one pixel of the object is saturated (or very close to) 

8 The object is truncated (too close to an image boundary) 

16 Objects aperture data are incomplete or corrupted 

32 Objects isophotal data are incomplete or corrupted 

64 A memory overflow occurred during deblending 

128 A memory overflow occurred during extraction 

this information was taken from the SExtractor documentation. 

These flags are produced for each band and are included in the DES 

DM database products. 


APPENDIX A: EXPLORATORY DATA ANALYSIS 

Given below are some examples of the types of object removed 
during each step of the selection process. 

Al Flags 

The DES flags were used as detailed in step one of section 3. In fig¬ 
ure Al we show examples of objects that fail only the first part of 
our flag cuts. As can be seen most of these objects are junk rather 
than valid candidates. The central cutout shows an object that is 
undetected in that band which are exactly the objects that we want 
to keep in the g and r bands as they correspond to high redshift 
quasar candidates. To ensure that candidates were not removed ob¬ 
jects were only removed if they failed both of the flag cuts. The sec¬ 
ond cut was made in the weight flags and the combination of these 
removed the junk objects but kept the faint, undetected sources in. 

The internal flags that are produced by SExtractor are: 

1 The object has neighbours, bright and close enough to signifi¬ 
cantly bias the MAG AUTO photometry, or bad pixels (more than 
10% of the integrated area affected) 


A2 g and r Band Detections 

In steps 5 and 6 a catalogue based method for removing objects 
that are detected in the bluer bands is used. Eirst objects with a 
magnitude brighter than 23 were removed. It was found that this 
did not reliably remove all detected objects. Some objects were in 
a particularly deep part of the survey and still had a clean detection, 
shown in figure A2, despite being fainter than 23rd magnitude. To 
remove these an error cut was used, sources with ag/r >0.1 were 
kept. 

The cut of ag/r >0.1 was used as it was empirically found 
from the data that this was a good balance between keeping a few 
faint objects and not throwing away objects that were falsely de¬ 
tected in g or r. Positive noise at the site of the object can cause 
SExtractor to register a magnitude for the object when there isn’t a 
clean detection. This is shown in the cutouts in figure A3. 

As high redshift quasars are unlikely to be detected in the g or 
r bands it is important that objects with a false detection in g or r 
are kept. 

Another problem that was found was that some objects had a 
g /^psf = 99 (which can mean a non detection) but were actually 
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Figure A3. This object has no discernible flux in the g band but a catalogue 
entry of ^fpgf = 24.88 ± 0.15 and has no detection in the r band. Despite 
appearing similar in both bands the object has been treated differently by 
the pipeline in each band; leading to a detection in one and not in the other. 
This demonstrates how the selection technique for dropouts has to take into 
account the multiple ways of treating objects at the detection threshold. 


q r i z Y 

• lelol»l4 


Figure A4. Here an object which has saturated in the g and r bands is shown. 
While it is obviously present in all the bands its PSF magnitudes are 99. This 
objects has ^fpgf = 99.0 ± 1.0, rpgf = 99.0 ± 1.0, Zpgf = 99.0 ±1.0 and 
= 14.7326 ± 0.0001 


very bright, saturated objeets. An example of this is shown in figure 
A4. 

The object in figure A4 has the same r band magnitude as the 
object shown in figure A3. To distinguish between these two cases 
a rough version of forced photometry was used as described in step 
9 of section 3. As the difference in flux values for these two cate¬ 
gories of object is very dramatic and the r band magnitudes were 
not required a complicated method of accurate forced photometry 
was not needed. 

A3 Cosmic Ray Removal 

A cosmic ray hit in the z band can cause an object to pass the selec¬ 
tion criteria as it looks like a very extreme i - z dropout. An example 
cosmic ray detection is shown in figure A5, this shows the irregular 
shape of the detection. 

The first step in removing these objects, which made up most 
of our sample after step 7 in section 3, was to require the objects to 
have a Y band magnitude. As the Y band data is shallower and more 
variable than the other bands no error cut was used. This left in 
objects that were false detections in the Y band (such as the object 



Figure A5. Cosmic ray hits in the z band create objects that pass the selec¬ 
tion due to looking like extremely red i - z dropouts. These cutouts show 
the odd shape of the cosmic ray as well as demonstrating that they are not 
smoothed out by the PSF. 


in figure A5 which has Ypgf = 20.54 ± 1.0) and the candidate list 
was still swamped with cosmic rays. To remove these image based 
cosmic ray removal was used. This step worked off the differences 
between adjacent pixels as it was found that real objects had much 
shallower gradients between adjacent pixels than cosmic rays. 

A4 SQL Used 

Here we present the SQL code used to query the DBS database. 
This code returns 838 rows from a database with 40,129,963 unique 
sources as shown in Steps l-8b in Table 1. 

SELECT 

* 

FROM 

SVAl_COADD_GRIZY 

WHERE 

/* Flag criteria */ 

(FLAGS_G = 0 OR FLAGS_WEiGHT_G = 0 OR 
FLAGS_G > 3) AND 

(FLAGS_R = 0 OR FLAGS_WEiGHT_R = 0 OR 
FLAGS_R > 3) AND 

(FLAGS_i = 0 OR FLAGS_WEiGHT_i = 0 OR 
FLAGS_i > 3) AND 

(FLAGS_Z = 0 OR FLAGS_WEiGHT_Z = 0 OR 
FLAGS_Z > 3) AND 

(FLAGS_Y = 0 OR FLAGS_WEiGHT_Y = 0 OR 
FLAGS_Y > 3) AND 

/-k Z band limits */ 

MAG_PSF_Z <= 21.5 AND MAGERR_PSF_Z < 0.1 
AND 

/-k Morphological criterion */ 

MAG_PSF_Z - MAG_MODEL_Z < 0.145 AND 

/* i-Z colour criterion */ 

MAG_PSF_i - MAG_PSF_Z > 1.694 AND 

/-k G and R dropout criteria */ 

MAG_PSF_G >23 AND MAG_PSF_R >23 AND 
MAGERR_PSF_G >0.1 AND 
MAGERR_PSF_R >0.1 AND 

/-k Z-Y colour criterion */ 

MAG_PSF_Z - MAG_PSF_Y < 0.5 AND 

/-k Y band detection criterion */ 

MAG_PSF_Y <> 99 AND 
MAG_PSF_Y <23 
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